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A scheme suggested in the literature to determine the symmetry-imposed shape of linear response 
tensors is revised and extended to allow for the treatment of more complex situations. The exten¬ 
ded scheme is applied to discuss the shape of the spin conductivity tensor for all magnetic space 
groups. This allows in particular investigating the character of longitudinal as well as transverse 
spin transport for arbitrary crystal structure and magnetic order that give rise e.g. to the spin Hall, 
Nernst and the spin-dependent Seebeck effects. In addition we draw attention to a new longitudinal 
spin transport phenomenon occurring in certain nonmagnetic solids. 
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I. INTRODUCTION 

The shape of a linear response tensor is of central im¬ 
portance as it decides whether a physical phenomenon 
may occur and what anisotropy may be expected for a 
solid with given crystal symmetry and magnetic order. 
A prominent and common example for this is the an¬ 
omalous Hall effect in ferromagnetic solids, that is con¬ 
nected with the non-zero anti-symmetric contributions 
to the electrical conductivity tensor. Accordingly, sev¬ 
eral schemes were suggested in the past to predict the 
shape of linear response tensors on the basis of group- 
theoretical arguments (for a corresponding review see for 
example Ref. 1). Among the various schemes suggested 
that of Kleiner 2-1 seems to be most convincing as it is 
starting from the expression for linear response tensors 
as given by Kubo’s linear response formalism and as it 
uses only the behavior of the involved operators under 
the appropriate space and time transformations of the 
relevant magnetic space group. A further appealing fea¬ 
ture of Kleiner’s scheme is that it does not make use of 
Onsager’s relations but allows to derive them in a most 
general way. 


Kleiner’s scheme was originally derived having re¬ 
sponse quantities in mind that are connected with the 
perturbation as well as the response represented by 
the components of a vector operator. A more general 
starting point is adopted in this contribution to allow 
the treatment of more complex situations. As a first 
simple application the tensors representing the charge 
and heat transport in response to an electric field and 
thermal gradient are considered. As a more complex 
transport quantity the corresponding spin conductivity 
is considered for all magnetic space groups. Among 
other things this allows the discussion of the transverse 
spin transport as occurring for the spin Hall ■ and 
spin Nernst ’' effects. In particular it is demonstrated 
that these effects may be discussed without use of 
spin-projected conductivities''’ 11 . 

II. SYMMETRY OF RESPONSE FUNCTIONS 

Within Kubo’s linear response formalism, the change 
of the expectation value of an observable Bi due to a 
time-dependent perturbation Aj can be expressed by the 
corresponding response function : 


TpiAj ( w > H) = J dt e~ iuit J ^ dA Tr ( p(U)AjBi(t + ihX; H)) 


(1) 


Here p(H) = e £*d H )/Tr(e ^( H )) is the density oper¬ 
ator for the unperturbed system, the operators Bi and Aj 
in the Heisenberg picture are assumed to be the Cartesian 
components of a corresponding vector operator and H is 
an external magnetic field. 

Eq. (1) was used by Kleiner 2 as the starting point to 
investigate the symmetry of the tensors t that describe 
the charge and heat transport due to an electric field 


or thermal gradient. Kleiner’s scheme, however, is quite 
general and can be easily extended to deal with more 
complex situations. In the following, Kleiner’s scheme 
will be adopted to the case when the observable is rep¬ 
resented by an operator product of the form BiCj , again 
with the operators C;, Bj, and Ak being the Cartesian 
components of a vector operator. In this case the corres¬ 
ponding response function is obviously given by: 
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4fc (w,H) = J df e~' ut dATr (^p(H)AkBi(t + ihX; H)Cj(t + ihX; H)^ 


'( B z C j )A k 


where by using the parenthesis in the symbol T (B-c )A k 
it is made clear that it is the observable and not the per¬ 
turbation consists in a product of two operators. The 
shape of the response tensor r in Eq. (2), i.e. the occur¬ 
rence and degeneracy of non-zero elements, has to reflect 
the symmetry of the investigated solid. This shape can 
be found by considering the impact of a symmetry op¬ 
eration of the corresponding space group on Eq. (2), as 
this will lead to an equation connecting the elements of 
t or possibly of a complementary tensor t' (see below). 
Collecting the restrictions imposed by all symmetry op¬ 
erations the shape of r is obtained. In this context it is 
important to note that the magnetic structure of the sys¬ 
tem, if present, has to be considered. In this case, the set 
of symmetry operations contains unitary pure spatial (u), 
but also anti-unitary symmetry operations (a). 

The general transformation properties of the operators 
Xi = Ai, Bi or Ci in Eq. (2) under unitary ( u ) and anti¬ 
unitary symmetry operations (a) can be written as: 


( 2 ) 


where D ! V - l (u) and D (X) (a) are the Wigner D-matrices 
corresponding to the operator X and operation u or a, 
respectively. The group properties are reflected by the 
following relations: 


D(uu') = D(u)D(u') 
D(aa') = D(a)D(a')* 


( 5 ) 

( 6 ) 


uXiU 1 = XjD^\u)jj 
j 

aXid- 1 = '52XjD(*\a)j i , 


( 3 ) 

( 4 ) 


For all unitary operations u the expression under the 
trace in Eq. (2) can be reformulated by cyclic permuta¬ 
tion and by inserting the factor u~ l u = 1: 


Tr (e-^^A k Bi{t + iftA; H )C j (t + ihX ; H)) = Tr (u^ue-^^u^uA^uB^t + ihX, H)u~ 1 uCi(t + ifiA; H)) 
^ite -/3 ^4 H )u -1 ^ (uBi(t + ihX,U)u~A (uCj(t + ihX, H)u -1 ^ . (7) 


= Tr 


The four expressions grouped in parenthesis can now be 
dealt with separately. The term containing Ak can be 
rewritten using Eq. (3). For the term containing Bj one 
has accordingly : 


I- 

with H u the transformed field 

uIF(H)u _1 = H{ H„) 


(9) 


uBi(t + ihX,H)u 1 = B m (t + ihX, H„ 

m 

D^\u) mi , 


( 8 ) 


connected with the operation u. For the term contain¬ 
ing Cj(t + ifiA,H) an analogous expression is obtained. 
Inserting these relations into Eq. (7) one obtains: 


Tr (e-^WAkBiit + iHX, H )C j {t + ihX, H)) = ^ Tr fe~^ H ^A,B m (t + i KX, H u )C n (t + ihX, H u ) 

Imn ^ 

D^\u) lk D^\u) ml D^\u) nj ] . (10) 


This equation must hold for any operators Ak, Bj and 
Ci, i.e. also in the special case Ak = Bj = Ci = 1, leading 


to: 


Tr (e"^^) = Tr . (11) 
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Inserting the two last equations into Eq. (2) for the gen¬ 
eral transport coefficients, one obtains the transforma¬ 
tion behavior of r under a unitary symmetry operation 
u : 

T (B i c j )A k ( w ’ H) = ^ T (B m C n )Ai ( w > 

Imn 

(i i) lk D& (u) mi D& (u) nj (12) 

A similar procedure can be applied for anti-unitary op¬ 
erators a that contain the time reversal T, i.e. that can be 

I 


decomposed as a = vT with v a unitary operator describ¬ 
ing a pure spatial operation. For anti-unitary operators 
cyclic permutation under the trace does not hold, but 
one may use the relation: 

Tr(aa / ) = [Tr(a/a)]* . (13) 

This expression can be used to transform Eq. (2) in a 
similar way as done for Eq. (7) leading to: 


Tr (e-Pn^AkBiit + ffiA, H )C 3 (t + ihX, H)) = Tr (a- 1 a*-P A Wa- 1 aA k a- 1 aB i (t + ffiA, H )a~ 1 aC j (t + ifiA, H)) 

1 * 

^ae _,3i ^ H )a _1 ^ ^aAfca -1 ^ [aBi{t + ihX, H)a -1 ^ (aCj(t + ffiA, H)a -1 ^ . (14) 


Of the four expressions in parenthesis, the second one which follows directly from the fact that a contains the 
is directly given by Eq. (4), while the first one can be time reversal operation and inserting these expressions 
rewritten by introducing H a via the definition into Eq. (14) one arrives at: 

aH{H)a- 1 = H{ H.) . (15) 

Expressing the last two terms according to 

aBi(t + ffiA, H)a -1 = B m (—t + iHX, H) 

m 

D^\a) mi , (16) 

I 

Tr (e-PX^AkBAt + ihX, H )CAt + ihX, H)') = V Tr 


D^{u)* lk D^\ur mi D^\ur nj . (17) 


one arrives at an expression that is completely analogous 
to Eq. (10): 


' / j \ 


Imn 


Using the relation : 

Tr (e-^iB(r)C(r)] = Tr (e~ pA A{-T)BC) (18) 


e -0H ( n a ) Ai z m (_ t + iftA) H a )C n (-t + ihX, H a ) 


Tr ^ e -^ (H) A k Bi(t + ihX, H)Cj(t + ihX, H)) = ^ Tr ( e -^ (H “><5Ji?t,iJ(i + ifiA, H a ) 

D^\uy ik D^\u)* mi D^\ur nj ) 


Imn 


(19) 


where 


(7 = 1, thus: 


(Tr L)* =Tr(L t ) 


( 20 ) 


Tr (e-^W) 


= Tr(e-^ H “)) . 


with L being a linear operator was used. Again, this 
equation must also hold for the special case A = B = 


( 21 ) 
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Finally, inserting all these relations one obtains the 
transformation behavior for r as 

= T AJ(cl fiL) ^g) 

Imn 

D&(a); k D&(a)* mi D&(ar nj ,( 22 ) 

which is the counter part of Eq. (12), but for anti-unitary 
operators a. 

It is important to note that in general the tensors 
T (BiC )A k an d r A t (s)c t ) are different objects represent¬ 
ing different response functions which are only interre¬ 
lated by Eq. (22). Accordingly, the symbols r and r' 
will be used below to distinguish them. Obviously, the 
two tensors r and r 7 coincide only if all operators and 
their adjoined ones are the same, i.e. A * = Bi and so on. 

Eqs. (12) and (22) relate the elements of the tensor r 
with all the elements of r and t 7 , respectively. As men¬ 
tioned above, these relations impose for each symmetry 
operation restrictions on the shape of r that allow to 
decide which elements have to be zero and which are de¬ 
generate. However, to find the final shape of r it is not 
necessary to derive restrictions for all symmetry opera¬ 
tions of the relevant space group. Instead, it is sufficient 
to use only a generating set of symmetry operations. Fi¬ 
nally, as was stressed by Kleiner 2 , for the application of 
Eqs. (12) and (22) it is not necessary to know the explicit 
form of the operators A it Bj and C^, but only their be¬ 
havior under a symmetry operation expressed by Eqs. (3) 
and (4). 

III. APPLICATIONS 

A. Symmetry operations and magnetic Laue 
groups 

For a periodic solid, the corresponding unitary sym¬ 
metry operations u can be represented by the Seitz 
symbol 12 : 

u = {R\t}, (23) 

where R describes a (proper or improper) rotation and t 
describes a translation. The application of this symmetry 
operation on a three dimensional vector v is defined as 

iiv = Rv + t , (24) 

where R is the three dimensional matrix representation 
of the rotation R and t a three dimensional translation 
vector. For an anti-unitary symmetry operation a, the 
time reversal operation T has to be considered in addition 
to the spatial symmetry operations. It can be included 
in the Seitz symbol according to: 

a = {R 1 1} T . (25) 

The transformation properties of a vector v under a de¬ 
pend now on its behavior under space inversion and time 


reversal. A vector that reverses its orientation under 
space inversion is called a spatial vector (or polar vec¬ 
tor), if it stays unaltered it is called a pseudo-vector or 
axial vector. 

Generally, the transformation of a vector field v(r) un¬ 


der an arbitrary symmetry operation s is given accord¬ 
ingly by: 

sv(r) = ±R v(s _1 r) , (26) 

where the sign is determined by the behavior of v(r) 
under time reversal T that may by part of s. On the 
other hand, a pseudo-vector field v(r) transforms as: 

sv(r) = ±det(R) R v(s _1 r) . (27) 

An example for this is the magnetic field H. As H 
changes sign under time reversal, the minus sign in 
Eq. (27) applies. In particular one has 

IH = +H (28) 

TH= H (29) 


for the application of space inversion I and time reversal 
T. In the following, we will use in parallel the symbols 1 
and l 7 for I and T, respectively. 

Taking into account the time reversal operation, the 
full symmetry of a periodic solid is represented by its 
magnetic space group Q that combines all symmetry op¬ 
erations of the type given in Eqs. (23) and (25). Alto¬ 
gether there are 1651 magnetic space groups that fall into 
three categories 1:1 : 

(a) G contains the time reversal operation T as an ele¬ 
ment, 

(b) G does not contain T at all, neither as a separate 
element nor in a combination, 

(c) Q contains T only in combination with another sym¬ 
metry element. 

Only nonmagnetic solids possess one of the 230 space 
groups of category a), while magnetically ordered solids 
belong either to category b) or c). Category b) con¬ 
sists of 230 space groups, isomorphic to the nonmagnetic 
space groups, and category c) combines the remaining 
1191 space groups. 

As the crystallographic magnetic point group of a peri¬ 
odic solid accounts for the translational symmetry de¬ 
termined by its Bravais lattice, it is sufficient to con¬ 
sider only the corresponding point group operations in¬ 
stead of the elements of its magnetic space group when 
dealing with Eqs. (12) and (22). Under certain condi¬ 
tions (see below) it is possible to restrict the considera¬ 
tion further to the corresponding magnetic Laue group 
of a solid, that is generated by adding the inversion op¬ 
eration / to the crystallographic magnetic point group. 
This conventional definition deviates from the older one 
used by Kleiner that derives the Laue group from the 
corresponding crystallographic point group by removing 
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from each improper rotation i? = Pr I its improper part 
I. For this reason we list in Tables I - III all magnetic 
point groups of the three categories together with their 
corresponding magnetic Laue group. The symbol in par¬ 
entheses gives in addition the magnetic Laue group as 
used by Kleiner"’. 


magnetic point group 

IT, I', IT 

21', ml', 2/ml', 2'/m, 2/m' 

2221', mm21', m'mm, 
mmml', m'm'm' 

41', 41', 4 /m', 4 /ml', 4 ' jm! 
4221', 4mml', 42ml', 

4m21', 4 /m'mm, 4' /m'm'm, 

4/mmml', 4'/m'mm , 4/m'm'm 
31', 3', 31' 

3121', 31ml', 3'lm, 3Tm', 31ml' 
3211', 3mll', 3'ml, 3'm'l, 3mll' 
61', 61', 6 '/m, 6/m', 6/ml' 

6221', 6mml', 6m2l', 


magnetic Laue group 

IT (T) 

2/ml' (21') 

mmml' (2221') 
4/ml' (41') 


4/mmml' (42221') 
31' (3') 

31ml' (3'2) 

3rall' (3'2) 

6/ml' (61') 


62ml', 6/m'mm, 6 '/mm'm, 

6'/mmm', 6/mmml', 6/m'm'm' 6/mmml' (6221') 

231', m'3', m3l' m3l' (23') 

4321', 43ml', m'3'm, m'3'm', m3ml' m3ml' (43'2) 


magnetic point group 

1 , 1 

2, m, 2/m 
222, mm2, mmm 
4, 4, 4/m 


magnetic Laue group 
1 (1) 

2/m (2) 
mmm (222) 

4/m (4) 


422, 4mm, 42m, 4m2, 4/mmm 4/mmm (422) 
3, 3 3 (3) 

312, 31m, 3lm 3lm (32) 

321, 3ml, 3m 1 3ml (32) 

6, 6, 6/m 6/m (6) 

622, 6mm, 6m2, 62m, 6/mmm 6/mmm (622) 
23, m3 m3 (23) 

432, 43m, m3m m3m (432) 


Table II: Magnetic point groups of category b) and their 
corresponding magnetic Laue group. In parentheses the 
magnetic Laue group according to its old definition 
used by Kleiner is given (see text). Because equivalent 
magnetic point group and Laue group symbols have not 
been removed (see text) there are 37 and 12 instead of 
32 and 11, respectively, entries. 


Crystallography 1 . 

B. Thermoelectric Coefficients 


Table I: Magnetic point groups of category a) and their 
corresponding magnetic Laue group. In parentheses the 
magnetic Laue group according to its old definition 
used by Kleiner is given (see text). Because equivalent 
magnetic point group and Laue group symbols have not 
been removed (see text) there are 62 and 12 instead of 
53 and 11, respectively, entries. 

Deriving the symbols for the magnetic point and Laue 
groups from those for the magnetic space groups, one 
may be led in some cases to two symbols that differ only 
concerning the sequence of the second and third generat¬ 
ors (one of these might be a dummy 1). As this depends 
on the chosen coordinate system and because the shape 
of the response tensor may depend on this choice, both 
symbols are listed although being completely equivalent. 
This applies to 31ml' and 3mll' of category a), 31m and 
3ml of category b) and 4 '/mm'm and 4'/mmm', 31m' 
and 3m'1 as well as 6 '/m'm'm and 6'/m'mm' of cat¬ 
egory c). Furthermore, it should be noted that for the 
magnetic Laue groups 2/ml' of category a), 2/m of cat¬ 
egory b) and 2'/m' of category c), the coordinate system 
has been chosen according to cell choice 1 of space group 
2/m as documented in the International Tables of X-ray 


Within linear response theory, the induced electric cur¬ 
rent density j and the heat current density q are given 
by 2 



with e = |e| the elementary charge and the electrochem¬ 
ical potential p which is related to the chemical poten¬ 
tial fa c and the electric potential p via 

p = p c - \e\p . (31) 

As explicitly demonstrated by Kleiner - ’ as well as below, 
the coefficients L, tJ satisfy Onsager relations of the form 

L n (H) = L n (-H) (32) 

L 22 (H) = L 22 (-H) (33) 

L 12 (H) = L£(-H). (34) 

Identifying the operators A t and Bi with one of the com¬ 
ponents of the electric current density operator j and 
the heat current density operator q and setting C/ = 1 
Eqs. (12) and (22) reduce to the expressions given by 
Kleiner to investigate the symmetry properties of the 
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magnetic point group 

magnetic Laue group 

to 

to 

2 '/m' (2') 

2 , 2 , 2, m'm2', m'm'2, m! 

m'm m'm'm (2'2'2) 

4', 4', 4'/m 

4'/m (4') 

4'2'2, 4'm'm, 4'2'm, 


4'm'2, 4'/mm'm 

4 '/mm'm (4'22') 

4'22', 4'mm', 4'2m', 


4'm2', 4'/mmm' 

4 '/mmm' (4'22') 

42'2', 4m'm', 42'm', 


4m'2', 4/mm'm' 

4/mm'm' (42'2') 

312', 31m', 3 lm' 

31m' (32') 

32'1, 3m'1, 3mT 

3m'1 (32') 

6', 6', 6 '/m' 

6 '/m' (6') 

6'2'2, 6 'm'm, 6'm'2, 


6'2'm, 6'/m'm'm 

6 '/m'm'm (6'22') 

6'22', 6 'mm', 6'm2', 


6'2 vrj, 6'/m'mm' 

6'/m'mm' (6'22') 

62'2', 6 m'm', 6m'2', 


62 'm!, Q/mm'm' 

Q/mm'm' (62'2') 

4' 32', 4'3m', m3m' 

m3m' (4'32') 


Table III: Magnetic point groups of category c) and 
their corresponding magnetic Laue group. In 
parentheses the magnetic Laue group according to its 
old definition used by Kleiner J is given (see text). 

Because equivalent magnetic point group and Laue 
group symbols have not been removed (see text) there 
are 52 and 13 instead of 37 and 10, respectively, entries. 

thermoelectric coefficients L, ,. His derivation will be re- 
peated her in a modified way as we use the conventional 
definition for the Laue group and as the results will be 
used later on. 

Expressing the electric current density operator j = 
— \e\v as a product of the electronic charge — |e| and the 
velocity operator v one can see that j transforms as a 
vector that changes sign under time reversal T and space 
inversion I: 

1 ji = ~j% (35) 

T ji = -ji ■ (36) 

The same relations apply for the heat current density op¬ 
erator qi ' . The corresponding 3x3 matrix representa¬ 

tion for a unitary operator u = {R\t} and an anti-unitary 
operator a = {R\t}T to be used in Eqs. (12) and (22) is: 

D (j) (u) = D w) (it) = D(i?) (37) 

D (j) (a) = D® (a) = -D (R) . (38) 


Eqs. (12) and (22) (with C) = 1) can be brought into a 
more convenient form by replacing every D by D = R 1 
and H by H u -i or H a -i, respectively. Thus, Eq. (12) for 
unitary operators u simplifies to: 

Tz iA] (w,H(i?)) = E T B fe A>>H) D(R)ki D(R)ij (39) 

kl 

and Eq. (22) for anti-unitary operators a to: 

r §iAj (co,-U(R)) = 5>,M>’H) D(R)ki D(R) tj (40) 

kl 

where 

H {R)i = Y J R ^{PR)H 3 ■ (41) 

i 

Here we used the fact that the matrices D(R)ij are real 
and that H is a pseudo-vector. A further simplification 
can be achieved by splitting R in a proper rotation Pr 
and the space inversion /, if it is contained in R. Expli¬ 
citly, this means that R = Pr if R is a proper rotation 
and R = PrI if R is an improper rotation. For proper ro¬ 
tations one has det(-R) = +1 while for improper rotations 
det(i?) = —1 holds. Because the space inversion amounts 
to a simple multiplication with — 13, this splitting can be 
expressed by: 

D (R) = det (R) D (Pr) . (42) 

Since the matrix D(l?) appears twice in Eq. (39) 
and (40), the two factors det(R) compensate each other, 
regardless whether R is a proper or an improper rotation. 
Thus, the final equation for the unitary operators is: 

r SiAj (u, H) = E T B k A t ( w > H )D(P R ) ki D(P R ) lj (43) 

kl 

and for anti-unit ary operators: 

= J2 T A l B k (^)D(PRhiD(P R )i j • (44) 

kl 

This splitting of R enables one to consider the symmetry 
property of the thermogalvanic coefficients of a solid on 
the basis of its magnetic Laue group instead of its mag¬ 
netic point group. This applies whether the conventional 
definition of the Laue group (see section III A) is applied 
or that used by Kleiner J . In the latter case the removal 
of the ineffective inversion / happens already when con¬ 
structing the Laue group. In the former case, one may 
add improper rotations R = Pr I where again I is inef¬ 
fective and Pr an element of both groups. Working only 
with the magnetic Laue group has the obvious advant¬ 
age that less cases have to be considered (see Table I - 
Table III) as there are only 32 magnetic Laue groups, 
while there are 122 different crystallographic magnetic 
point groups. 

On the basis of Eqs. (43) and (44) it is now rather 
straightforward to give explicit forms for the response 
tensors L ZJ in Eq. (30). For this purpose the abbrevi¬ 
ations Tij = T Aits ., r-j = t & . a . and a i:j = t a . a . will be 
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used, where A and B can stand for j or q. Accordingly, 
t and r 7 represent either L 12 or L 21 or the other way 
around, and g_ represents L n or L 22 , respectively, that 
obviously have to have the same structure. It is interest¬ 
ing to note that Eq. (44) can lead to restrictions on the 
tensor elements in addition to those imposed by Eq. (43). 
These hold even for the tensors of type r 7 . 

In the case of a magnetically ordered solid having a 
magnetic space group of category b) the restrictions to 
the shape of the thermogalvanic tensors result only from 
the application of Eq. (43) as there are no anti-unitary 
operations. As a consequence, all tensors <x, r and r 7 
have the same shape. Accordingly, only the shape of r is 
given in Table V, that is in full agreement with Kleiner’s 
Table IV 2 . 

For magnetic space groups belonging to category a) 
or category c) Eq. (44) has to be applied in addition to 
Eq. (43). In general, this leads to different symmetry re¬ 
strictions for the tensors of type r 7 and g_. The resulting 
shape of the tensors for category a) is given in Table 
IV. These results agree with those given by Kleiner’s 
Table V 2 , apart from those for the Laue groups 31', 
4/ml 7 and 6/ml'. Since the magnetic Laue groups in 
category a) differ from those in b) only by the time- 
reversal l 7 as an element of its own, the tensor shapes 
in Table IV alternatively can be deduced from those 
in Table V simply by considering in addition the effect 
of V. In case of <x this can lead to additional restric¬ 
tions (degeneracies and zero elements) since <x 7 = <r. 
For the thermoelectric tensor on the other hand, this 
just states the usual Onsager relations as expressed by 
tL(H) = Tji{— H) (see Eq. (34)). Table VI gives the res¬ 
ults for category c) that are in full agreement with those 
given by Kleiner’s Table VI 2 . Obviously, the results 


magnetic 
Laue group 


IT 


2/ml 7 

mmm l 7 

31', 4/ml 7 , 
6/ml 7 

31ml 7 , 3mll 7 , 
4/mmml 7 , 
6/mmml' 

m3l 7 , m3ml 7 


0 


Txx Tyx Tzx 
Txy Tyy Tzy 
Txz Tyz Tzz 
Txx 0 T Z x 
0 Tyy 0 
Txz 0 T zz 
Txx 0 0 

0 Tyy 

0 0 T; 

Txx Txy 0 
Txy Txx 0 
0 0 T z - 

Txx 0 0 

0 Txx 0 
0 0 T zz 

Txx 0 0 

0 Txx 0 
0 0 Txx 


(Txx (Txy (Txz 
(Txy (Tyy (Tyz 
(Txz Tyz (Tzz 
(Txx 0 (Txz 

0 (Tyy 0 



0 

0 

Vyy 

0 

0 

(Txx 

0 

0 

(Txx 

0 

0 

(Txx 

0 



Table IV: Tensor forms for magnetic Laue groups of 
category a). 


presented in Tables IV - VI fulfill the Onsager relations 


magnetic Laue group 

I 

2/m 

mmm 

3, 4/m, 6/m 

3lm, 3ml, 4/mmm, 6/mmm 
m3, m3m 


xx I xy 1 xz 
r yx Tyy Tyz 


1 zx • zy • zz 


Txx 0 T X z 
0 Tyy o 
T zx 0 Tzz 
T X x 0 0 

0 Tyy 0 
0 0 Tzz 

Txx Txy 0 

~Txy Txx 0 

0 0 Tzz 

Txx 0 0 

0 Txx 0 
0 0 Tz ; 

Txx 0 0 

0 Txx 0 
0 0 Txx 


Table V: Tensor forms for magnetic Laue groups of 
category b). 


given by Eqs. (32) to (34) that are not postulated a pri¬ 
ori. 

Kleiner’s scheme was applied here to derive the shape 
of the tensors representing homogeneous bulk systems. 
However, it may also be applied to investigate the sym¬ 
metry restrictions on the so-called layer-resolved con¬ 
ductivity tensor a_ IJ with I and J labeling atomic lay¬ 
ers of a two-dimensional periodic system . This concept 
may be used for example in the context of electrical trans¬ 
port in layered GMR systems ,1, or magneto-optical 
properties of surface systems 1 2 . Another extension of 
Kleiner’s scheme is the discussion of non-linear effects 11 '. 


C. Shape of the spin conductivity tensor 

Spin transport as reflected for example by the spin Hall 
effect is usually described by use of the spin conductiv¬ 
ity tensor that gives the current density along direc¬ 
tion i for the spin polarization with respect to the fc-axis 
induced by an electrical field along the j-axis. Within 
a single-particle description of the electronic structure 
the Kubo-formalism leads for cT to an expression ana¬ 
logous to the Kubo-Bastin equation ’ 1 for the electrical 
conductivity 22 ’ 23 : 

4 = w/ dsm 
Tr 

~Ji S(E - = - (45) 

Here H is the Hamiltonian of the system, G + (E) and 
G~(E) are the corresponding retarded and advanced 
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~Txy Txx 0 
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0 Txx 0 
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0 T x x 0 
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Txx 0 0 

0 T xx 0 
0 0 Txx 


Table VI: Tensor forms for magnetic Laue groups of 
category c). The tensor forms for the groups A'/rnrn'rn 
and 4' /mmm' are related to each other by a rotation of 
the coordinate system around the principal (z) axis by 

7t/4. 


Green functions and jj is the ordinary electrical cur¬ 
rent density operator. A straightforward definition for 
the spin current density operator J// = \{vi,cr k } con¬ 
sists in the anti-commutator of the conventional velocity 
operator and the Pauli spin matrix a k - 21 As the spin 
conductivity is caused by spin-orbit coupling a coher¬ 
ent relativistic implementation of Eq. (45) seems to be 
more appropriate. This implies that the electrical cur¬ 
rent density operator jj = — e|cay is expressed in terms 
of the 4x4 Dirac a-matrices 25 . A corresponding ex¬ 
pression for the spin current density operator = Tkji 
was suggested by Vernes et al. 2h that involves the spatial 
part Ik of the spin polarization operator introduced by 
Bargmann and Wigner 

ti = PVi~ ^7slIi • (46) 

Here /?, 75, X, are again standard 4x4 Dirac-matrices, 
m is the electron mass and Hi stands for the kinetic 
momentum 2 . In fact this approach was adopted by Low- 
itzer et al . I2,2:! when dealing with the spin Hall effect of 
disordered alloys. However, as mentioned above, for an 
investigation of the shape of a response tensor the expli¬ 
cit expressions for the involved operators are not relev¬ 
ant but only their behavior under symmetry operations. 


Both definitions of J 1 / given above, consist of a combina¬ 
tion of the velocity operator jj with an operator that rep¬ 
resents the spin polarization of an electron. In contrast 
to ji (see Eq. (38)), the latter one (e.g. 7)) transforms as 
a pseudo-vector which changes sign under time reversal. 
Accordingly, one has for the transformation matrices 

D (Ti) (m) = det(ii)D(J?) (47) 

D (Ti) (a) = - det(i?)D(fl) , (48) 

corresponding to Eqs.(5) and (6). 

Identifying now Aj = j), Bj = ji and C k = T k in 
Eqs. (12) and (22) one finds the behavior of oA un der 
unitary transformations 

4 = E det W D ( R U °lm (49) 

Imn 

and under anti-unitary transformations 
4 = - ^imn det (R) D{R)u D{R) m j D(R) nk a'Z , (50) 

respectively. In analogy to the treatment of thermoelec¬ 
tric coefficients presented above one may again split the 
rotation R into its proper part P R and, if present, im¬ 
proper part as given in Eq. (42). The resulting equation 
for unitary transformations is then: 

4 = E^( p ^ D ( p R)mj D{P R ) nk det (A) 4 <T? m (51) 

Imn 

= Y, D(P r ) u D{P R ) mj D(P R ) nk af m (52) 

Imn 

and 

4 = - E D ( P rU D{P R ) nk a'Z (53) 

Imn 

for anti-unitary transformations, respectively. As a con¬ 
sequence, as found for the thermoelectric coefficients by 
Kleiner 2 also for the spin conductivity tensor it is suffi¬ 
cient to consider the magnetic Laue group of the solid. 

Using Eqs. (52) and (53) the shape of the spin con¬ 
ductivity tensor was determined with the results given 
in the left column of Tables VII - IX for magnetic Laue 
group of category a) - c). 

It should be noted that these constitute the equivalent 
to the generalized Onsager relations derived by Kleiner 
for the t' tensors. 

Because j and q have the same transformation prop¬ 
erties and because the tensors 7 ^ c )A k an d t a\b‘ [ c + ) 
in Eq. (22) are different objects in both cases, the tensor 
shapes for tensors describing the connection between heat 
currents and spin currents have exactly the same shape 
as those tabulated in Tab. VII, VIII and IX. 

For convenience, it is possible to alter the notation of 
these symmetry-restricted matrices in such a way that 
the symmetry of the tensor is easier to recognize at first 
sight. However, this reduction leads to the loss of the 
specific meaning, i.e. the generalized Onsager relations, 
contained in the tensors cr' k . The reduced tensors are 
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tabulated for category a), b) and c) in the right column 
of Tables VII, VIII and IX, respectively. 

As discussed in the context of the charge and heat 
current in response to an electric field the corresponding 
operators ji and qi have the same symmetry properties. 


As a consequence the tensors L n and L 22 in Eq. (30) 
have the same shape given by rr in Tables IV - VI. For 
the same reason the tensor representing the spin current 
induced by an thermal gradient has the same shape as 
that connected with an electric field with both given by 
Tables VII - IX. 
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0 0 0 
0 0 -al z 
) 

_x y x 

GXX V XX Qxz 

n-y —rr x —rr y 

u XX u XX u xz 
0 x z X -0 L 0 

ol* o 0 

_x _y 
(T —(T y 
^ XX ^ X 

0 -0zx 0 

0 0 XX 0 

*lx 0 -0 y 

0 -0 y zx 0 

0 0 0 
0 0 0 X i 

O0 y z o* 

0 0 0 
0 0 0 XJ 
0 0 xy 0 


cr y cr y cr y 

^ xx ^ xy ^ xz 

g v a y a y 

yx yy w yz 

a y a y „y 

zx zy zz 


0 0yy 0 

r y xx 0 0 « 

o 0 0 XZ 
0 0 0 
0ZX 0 0 
0 0 0 y z 
0 0 0 X Z 
T zx 0zx 0 
0 0 0 y z 
0 0 0 


' XX XX '-'xz 

x _ y x 

' XX "XX V xz 


' zx ^zx 


0 0 


0 ol 


0 -ol* 0 
0 y 0 0 


0 0 0 
0 x y 0 0 

0 O-0 X 
0 0 0 
0 X y 0 0 


0 0 xy 0 

z rv z 
'-'yx '-'yz 

0 0 Z zy 0 

0 0 X y O' 

0yx 0 0 

0 0 0 
2 z rv 
^xx '-'xy ^ 

_ 2 2 rv 

Vxy v xx ^ 


0 0 o. 

r£y( 
0 0 
0. 
0 
( 

0 0 a; 

r£y( 
0 0 


0 0 X y O' 

_ 2 

v xy 

0 0 0 


_ 2 2 rv 

xy v xx ^ 


0 0 X y O' 


0 xy 

0 0 0 . 

0 0%y 0 

— 0 x y 0 0 

0 0 0 
0 0 X y 0 
0 V XZ 0 0 

0 0 0 
0 0%y O' 

-0xy 0 0 

0 0 0 


Table VII: Polarization tensor forms and reduced polarization tensor forms for magnetic Laue groups of category a). 
The tensor forms for the groups 31ml' and 3?nll' are related to each other by a rotation of the coordinate system 

around the principal (z) axis by 7 t/2. 
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Table VIII: Polarization tensor forms and (identical) reduced polarization tensor forms for magnetic Laue groups of 
category b). The tensor forms for the groups 31m- and 3ml are related to each other by a rotation of the coordinate 

system around the principal (z) axis by n/2. 



11 


magnetic 

Laue 

group 

2 '/m' 


4'/m 
4' /mm'm 

4 ' /mmm! 

A/mm'm ', 
6/mm'm' 

31 m' 

3m'1 
6'/m' 


6'/m'm'm 


6'/m'm-m' 
m3m' 


— —(T 

rr rr n-.'ii rr 


0 0 a* 

o o -al. 

z _ z pi 

' xx ®xy ^ 

0 0 —al z 

0 0 a* yz 

Z Z rv 

. rr..,u 


0 0 G% z 

0 0 -al 

Z _ Z ri 

a :rx a x y u 

0 0 0 

0 0 a% x 

OCyx 0 

0 0 a% z 

0 0 -al 

g.z _ q 


XX 

— <J X G y G ’ 

^ xx ^ xx ^ xz 

rrV rr x — rr y 
^ XX XX ^ xz 

Z _ Z pi 

77 XX & X11 


J xx ^ xy 

_x _y _x 

rr — ( t a rr 

^ xx ^ xx ^ xz 
— /jV — (j x —cjV 

^ xx ^ xx ^ xz 

Z Z pv 

77 x.x. 77 X y 


J XX 
X 


' XX ^ XX 


(T —G 

XX X 


(J a —G — G a 

z Z pv 

77 xx 77 X y 'J 

: T 0 0 


0 -a%y 0 


0 ctL 0 

rV 


XX 

0 al 


0 —a% y 0 


0 0 0 
0 0 a 
n^y 


xy 


^ yx 

—a* 


G yy G yz 

-rr y (T y 
u yy u yz 


yx 

0 

0 

_ z 

'-'yx 

o 

0 

z 

77 xy 


0 ~a* z 
0 a y yz 

Oyy 0 

o al z 

0_ -of, 
ala-. 0 


0 of, 
0 -of 

-of, o 

I 0 of, 

I 0 0 


a% v 0 0 


0 
0 

(7 rpq. 7. 


xy 

V 

XX 


(7, 

(7', 

<7\ 

—a. 

—g. 

z 

xy 


’ XX 

.z 

xy 

y 

XX 


v ^ xz 

o of, 

XX 0 

y 

XX xz 

y 

xx xz 

0 


77 XX 77. 

— (j y ^, (7 

a. 

— g~ a: 


’ XX 


XX 

z 

77x11 77. 


XX ^ xz 

0 


0 

77 x 


— rr rr y 

xx X 

0 0 

0 0 

0 G y r 


0 — c 

T z 

’xy 

0 

0 


0 0 
o of 


0 0 
0 0 


— G G —G 

7.x. 7.1! 7. 


„x ^.x r\ 

77ZX 77zy U 

-a L of„ 0 
0 0 o* 


— G 


v a y 

zy zx 
_x x 

77 Z n ~77 


0 

zy vzx 0 

0 0 -al 


0 0 


a ly 0 0 
0 0 0 
x y 
<Tzx a zx 


0 


^ zx zx 

0 0 a 


^ zx zx ^ 

0 0 a z zl 

o X zx of, 0 
-of, of, 0 
0 0 al. 


/ (jV 

/ ^ XX 

X 

77 xx 

rrlz \ 

f-rrlx 

G y 
u zx 

0 \ 

1 rr x 

1 77 xx 

— G y ~ 

^ XX 

-rrlz ) 

-rr 'lx ~ 

77 zx 

o. 

\ 77 X y 

z 

77 xx 

0 J 

\ o 

0 

-rr zz ) 


0 o“ T 0 
-of, C 
0 0 0 . 
0 of, 0 

-of, 6 o 
0 0 0 
o a% v 0 

rSz 0 0 
0 0 0 


: a^ y a^ z 

: a zy a z z 


0 of 


0 a: 


0 of 


0 

0 

X 

G. 

o 
0 

^•-of. 0 

0 0 0 
0 0 a, 

0 


0 a%z 
x n 

'zx 77 zy U 

0 of, 
0 -of 


y z 

0 


' zy 

0 a. 
0 -c 

-of, ( 

G y G 

— G^.-,. — G. 

— & 

G. 


xz 

y 

XX ^ xz 

7 ZX 0 


xz 

y 

XX ^ xz 


— G —G 

x.x. X. 


tV 


— G —G 
xx KJ x 

-of, 0 
0 0 
0 -of,-of 
0 -of, 

0 


0 

of, 0 
0 —of 
-of, 0 
I 0 0 

I 0 aly 
I of. 0 


a v a 
^ yx ^' 

Vof, a 


zy 

( 0 0 of 2 ^ 

0 0 of 2 

Vof, ofy 0 ) 

( 0 0 al z \ 

0 0 a% z 

Vof, of y 0 ) 
0 0 al z 

0 0 -of 

of,-of, 0 
/ 0 0of 2 \ 

0 0 0 
\of, 0 0 ) 


XZ 


( 0 0 a. 

0 0 a: 

\ -V —X 

\Gzx Gzx 

G y 
^ xx 


X 


^ XX 

— V _u; 

&ZX &ZX 

-erf 

G, 


-lz\ 

xz 

0 J 

^.x ^.y 

- G G y 
xx ^ xz 

G% z 

o' 


G y 
^ xx 


G y 
^ xx 


G xx 

G y 
KJ zx 


G y 
u xz 


XX 

■ y G X 

XX ^ xz 
X 

°zx 0 


y • 

rj y —rr 

xx : 


Jl, 

v xx 

y a 


XX 

G y 
^ xx 


0 

_ x 

^ xx 


G y 
u xz 

G% z 

0 

G y 
u xz 


G V 
u zx 

{G V 

' ^ XX 


G V > 
u xz 


zx 

_rr x 

Oxx 

0 
0 
0 

0 -of, 0 
of, 0 0 

( 0 Oof z \ 
0 0 0 
Vof^O 0 J 


^7^x G i^y G T£ z 
& zx & zy & zz 

z z p| 

^7xx 77xy ^ 

Z Z pj 
7* yx T' yy ^ 

0 0 of, 

rr xx rr xy 0 

Z Z Q 

Oy, a yy u 

o o oL 


rr xy o,, 0 
0 0 0 


0 aly O' 
rr z yx 0 0 
0 0 0 

Z Z n 

rr xx a x y n 

_ Z Z rv 

rr xy a xx ri 
0 0 of, 

z z ri 

o,, a xy U 

_ z z q 

rr xy o,, u 
0 0 of; 

z z Q 

a xx a xy U 

_ z z Q 

rr xy rr xx ^ 
0 0 Of; 

Z Z Q 

a xx r) xy r) 
_ z z r» 

rr xy o,, v 


0 0 o, 

r Z xy ( 
0 0 


0 Of y O' 


rr xy 
0 0 0 . 

0 of „ 0 

-ofy 0 0 

0 0 0 

o ofyO 
rfz 
0 0 0 


0 0 


Table IX: Polarization tensor forms and reduced polarization tensor forms for magnetic Laue groups of category c). 
Note that the reduced tensor forms for the groups m'm'm and 4 '/m as well as for the groups 31TO', 3m'l, and 6'/m' 
are identical. Moreover those of 7/mm'rn and 4 '/mmm', of 31?n/ and 3m'\ as well as of 6'/m'm'm and 6' / m!mm' 
are (pairwise) related to each other by a rotation of the coordinate system around the principal (z) axis by 7r/4, tv/2, 

and 7 t/ 2, respectively. 


Obviously, the occurrence of anti-symmetric off- 
diagonal elements in the tensor <r s (s=x, y, z) in Table 
VII implies that the transverse spin Hall effect is, in prin¬ 
ciple, allowed by symmetry in any paramagnetic solid. 
However, one has to stress that in case of the magnetic 
Laue groups 1, 2/m, and mmml' the shape of the tensor 
is not purely anti-symmetric. The same is true for a fer¬ 
romagnetic solid according to Tables VIII and IX, i.e. 
the spin Hall and Nernst effects are symmetry-allowed 
in any magnetic solid as well (again not all cases show 
purely anti-symmetric elements). Considering as an ex¬ 
ample a ferromagnetic cubic solid with the magnetic Laue 
group 4/mm'm' (e.g. bcc-Fe or fcc-Ni with the magnet¬ 
ization along z-direction) its spin conductivity tensor is 


very different from the form of its nonmagnetic counter¬ 
part with m3ml'. For the nonmagnetic case only the ele¬ 
ments tr*) with i zjz j zfz k zjz i are non-zero. In addition, 
these are the same for a cyclic permutation of (i,j, k) and 
change the sign for an anti-cyclic one. For the ferromag¬ 
netic case additional off-diagonal elements may appear, 
with the degeneracies depending on the spin projection 
component fc, and the tensors are no longer purely anti¬ 
symmetric. In particular one notes that there are diag¬ 
onal elements that imply the occurrence of a longitudinal 
spin current induced by an electric field that in general 
will depend on whether the electric field is along (erf. ) 
or perpendicular (. a z y = a yy ) the magnetization. These 
tensor elements are obviously responsible for the occur- 
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rence of the spin-dependent Seebeck effect' . Interest¬ 
ingly, for a nonmagnetic solid there are several magnetic 
space groups that imply a non-vanishing diagonal tensor 
element oA, i.e. a longitudinal current along the direction 
of the applied electric field or thermal gradient. This was 
demonstrated recently by corresponding numerical work 
on nonmagnetic (Aui_ x Pt x ) 4 Sc showing that the longit¬ 
udinal spin conductivity can be comparable in magnitude 
to the transverse spin Hall conductivity. ' 1 

D. Implementation 

The symmetry restrictions imposed on the thermogal- 
vanic tensors by Eqs. (43) and (44) as well as on the 
spin conductivity tensor by Eqs. (52) and (53), respect¬ 
ively, were determined by means of a Python script that 
is based on the Computational Crystallography Toolbox , 
cctbx . Although this library provides support only for 
the nonmagnetic crystallographic operations, it is also of 
great value when dealing with magnetic solids. To de¬ 
termine the magnetic space group of a solid all possible 
magnetic space groups are simply scanned through and 
checked which fits to the system under investigation. The 
corresponding symmetry operations are taken from the 
magneticspacegroupdatafilemagnetic_data.txt 
Once the magnetic point group has been found, the u 
and a operators needed for an application of Eqs. (43) 
and (44) or Eqs. (52) and (53), respectively, are fixed. 
Going through all elements of the magnetic point group 
leads to a set of connecting equations between the tensor 
elements which can then be solved to get the shape of 
the tensor. For these symbolic calculations the SymPy 
library is used. Although in principle the generators of 
a magnetic point group are sufficient to obtain all sym¬ 
metry restrictions, it turned out to be more convenient 
to apply all symmetry operations since the cctbx library 
and the magnetic space group tables do not provide a set 
of generators. 

Finally, it should be mentioned that the results for 
the spin conductivity tensor <x 2 for the spin polarization 
along the 0 -axis have been checked against the output 
of the SPRKKR program package"" that allows to cal¬ 


culate this tensor on the basis of the relativistic Kubo 
formalism . For all investigated magnetic Laue groups 
of category a) (11', mml', 2/ml', A/ml', 4/mmml', 
6/mmml', m3ml'), b) (4/m), and c) (2 '/m!, m'm'm , 
A/mm'm', 3m' 1, 6/mm'm') the numerical results for a_ z 
were found to be completely in line with the analytical 
predictions given in Tables VII - IX. 

IV. SUMMARY 

Kleiner’s scheme to determine the shape of a linear re¬ 
sponse tensor have been extended to deal with more com¬ 
plex situations. The scheme has been used to revise the 
shape of the electric charge and heat conductivity tensors 
for all magnetic space groups. It was demonstrated that 
for this only the magnetic Laue group of a solid is rel¬ 
evant. This also holds for the spin conductivity tensor, 
that is used among other to discuss the longitudinal spin- 
dependent Seebeck effect as well as the transverse spin 
Hall and Nernst effects. Results for all magnetic space 
groups are presented in an easily accessible way, by giv¬ 
ing in addition to the tensors g/ k containing the general¬ 
ized Onsager relations also the reduced tensor forms <x fc . 
Furhermore, the axis conventions of the space groups are 
preserved when constructing the magnetic Laue groups 
and therefore, although redundant, the tensor forms are 
given in both coordinate systems whenever there is an 
ambiguity. Interestingly, several magnetic Laue groups 
for nonmagnetic solids were identified that should show 
a new longitudinal spin transport phenomenon. 9 Finally, 
it should be stressed that the scheme presented here can 
be applied straightforwardly to any other response func¬ 
tion. Examples relevant for spintronics and related fields 
are the response tensors representing spin-orbit torque, 
Gilbert damping or the Edelstein effect." 
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